RNA three-dimensional structure drives the sequence organization of potato spindle tuber viroid quasispecies.

RNA viruses and viroids exist and evolve as quasispecies due to error-prone replication. Quasispecies consist of a few dominant master sequences alongside numerous variants that contribute to genetic diversity. Upon environmental changes, certain variants within quasispecies have the potential to become the dominant sequences, leading to the emergence of novel infectious strains. However, the emergence of new infectious variants remains unpredictable. Using mutant pools prepared by saturation mutagenesis of selected stem and loop regions, our study of potato spindle tuber viroid (PSTVd) demonstrates that mutants forming local three-dimensional (3D) structures similar to the wild type (WT) are more likely to accumulate in PSTVd quasispecies. The selection mechanisms underlying this biased accumulation are likely associated with cell-to-cell movement and long-distance trafficking. Moreover, certain trafficking-defective PSTVd mutants can be spread by functional sister genomes in the quasispecies. Our study reveals that the RNA 3D structure of stems and loops constrains the evolution of viroid quasispecies. Mutants with a structure similar to WT have a higher likelihood of being maintained within the quasispecies and can potentially give rise to novel infectious variants. These findings emphasize the potential of targeting RNA 3D structure as a more robust approach to defend against viroid infections.

Immune System Deficiencies Do Not Alter SARS-CoV-2 Evolutionary Rate but Favour the Emergence of Mutations by Extending Viral Persistence.

During the COVID-19 pandemic, immunosuppressed patients showed prolonged SARS-CoV-2 infections, with several studies reporting the accumulation of mutations in the viral genome. The weakened immune system present in these individuals, along with the effect of antiviral therapies, are thought to create a favourable environment for intra-host viral evolution and have been linked to the emergence of new viral variants which strongly challenged containment measures and some therapeutic treatments. To assess whether impaired immunity could lead to the increased instability of viral genomes, longitudinal nasopharyngeal swabs were collected from eight immunocompromised patients and fourteen non-immunocompromised subjects, all undergoing SARS-CoV-2 infection. Intra-host viral evolution was compared between the two groups through deep sequencing, exploiting a probe-based enrichment method to minimise the possibility of artefactual mutations commonly generated in amplicon-based methods, which heavily rely on PCR amplification. Although, as expected, immunocompromised patients experienced significantly longer infections, the acquisition of novel intra-host viral mutations was similar between the two groups. Moreover, a thorough analysis of viral quasispecies showed that the variability of viral populations in the two groups is comparable not only at the consensus level, but also when considering low-frequency mutations. This study suggests that a compromised immune system alone does not affect SARS-CoV-2 within-host genomic variability.

Quasispecies productivity.

The quasispecies theory is a helpful concept in the explanation of RNA virus evolution and behaviour, with a relevant impact on methods used to fight viral diseases. It has undergone some adaptations to integrate new empirical data, especially the non-deterministic nature of mutagenesis, and the variety of behavioural motifs in cooperation, competition, communication, innovation, integration, and exaptation. Also, the consortial structure of quasispecies with complementary roles of memory genomes of minority populations better fits the empirical data than did the original concept of a master sequence and its mutant spectra. The high productivity of quasispecies variants generates unique sequences that never existed before and will never exist again. In the present essay, we underline that such sequences represent really new ontological entities, not just error copies of previous ones. Their primary unique property, the incredible variant production, is suggested here as quasispecies productivity, which replaces the error-replication narrative to better fit into a new relationship between mankind and living nature in the twenty-first century.

The dynamic variation position and predominant quasispecies of hepatitis B virus: Novel predictors of early hepatocarcinoma.

To find the predictors of early HCC based on the dynamic changes of HBV quasispecies, this study utilizing the second-generation sequencing (NGS) and high-order multiplex droplet digital PCR (ddPCR) technology to examine the HBV quasispecies in serum of total 247 subjects recruited from high-incidence area of HCC. In the discovery stage, 15 non-synonymous Single Nucleotide Polymorphisms (SNPs) with higher variant proportion in HCC case group were founded (all P<0.05). Furthermore, the variant proportions in some of these SNPs were observed changing regularly within 5 years before the onset of HCC, and 5 of them located in HBX, 2 in HBS and 2 in HBC. The HBV predominant quasispecies and their consensus sequences were identified by genetic evolution analysis, in which the high HBS and HBC quasispecies heterogeneity were found associated with the forming of multifarious quasispecies clones, and the HBX gene had the highest proportion of predominant quasispecies (46.7 % in HBX vs 12.7 % and 13.8 % in HBS and HBC respectively) with the key variations (G1512A, A1630G, T1753C/G/A, A1762T and G1764A) determined. In the validation stage, we confirmed that the combined double mutations of G1512A+A1630G, A1762T+G1764A, and the combined triple mutations of T1753C/G/A + A1762T+G1764A, all expressed higher in early HCC cases when comparing with control group (all P<0.05). We also demonstrated the advantages of ddPCR using in multi-variations detection in large-sample for early HCC surveillance and screening. So we think that the dynamic of key HBV variation positions and their different combinations determined by quasispecies anlysis in this study can act as the novel predictors of early hepatocarcinoma and suitable to popularize and apply in HCC screening.

Cell-cell communication and initial population composition shape the structure of potato spindle tuber viroid quasispecies.

RNA viruses and viroids replicate with high mutation rates, forming quasispecies, population of variants centered around dominant sequences. The mechanisms governing quasispecies remain unclear. Plasmodesmata regulate viroid movement and were hypothesized to impact viroid quasispecies. Here, we sequenced the progeny of potato spindle tuber viroid intermediate (PSTVd-I) strain from mature guard cells lacking plasmodesmal connections and from in vitro-cultivated mesophyll cell protoplasts from systemic leaves of early-infected tomato (Solanum lycopersicum) plants. Remarkably, more variants accumulated in guard cells compared to whole leaves. Similarly, after extended cell culture, we observed more variants in cultivated mesophyll protoplasts. Coinfection and single-cell sequencing experiments demonstrated that the same plant cell can be infected multiple times by the same or different PSTVd sequences. To study the impact of initial population composition on PSTVd-I quasispecies, we conducted coinfections with PSTVd-I and variants. Two inoculum ratios (10:1 or 1:10) established quasispecies with or without PSTVd-I as the master sequence. In the absence of the master sequence, the percentage of novel variants initially increased. Moreover, a 1:1 PSTVd-I/variant RNA ratio resulted in PSTVd-I dominating (>50%), while the variants reached 20%. After PSTVd-I-only infection, the variants reached around 10%, while after variant-only infection, the variants were significantly more than 10%. These results emphasize the role of cell-to-cell communication and initial population composition in shaping PSTVd quasispecies.

Modelling the structure and evolution of cultural information as quasispecies.

We present a novel mode of cultural evolution whereby some forms of transmission may be modelled as quasispecies. The model incorporates the effect of high rates of error in certain forms of communication; while also building on the structural similarities between biological molecules and written language. Firstly, both written language and key biological molecules, such as RNA and proteins, are modular. Within these molecules, structural domains may be recombined, while retaining their function. Likewise, sentences are structured as combinations of clauses, in which each clause contains a domain of information. The clausal structure permits the recombination of information to adopt different meanings, while allowing each unit to retain its identity. Secondly, by virtue of intrinsically-high error rates, we show that some, but not all, aspects of communicated culture information exists as rapidly evolving clouds within the population. These clouds of cultural information behave as quasispecies, which we model with varying mutation rates and suitable selection coefficients. We then integrate these ideas with the application of Shannon Diversity Index to produce a more holistic view of culture that is centred on the evolution of its information. Re-imagining culture, as evolving clouds of information, unifies the mode in which information is stored culturally and biologically, and opens up new avenues of comparative analysis.

Puzzles, challenges, and information reservoir of SARS-CoV-2 quasispecies.

Upon the emergence of SARS-CoV-2 in the human population, it was conjectured that for this coronavirus the dynamic intra-host heterogeneity typical of RNA viruses would be toned down. Nothing of this sort is observed. Here we review the main observations on the complexity and diverse composition of SARS-CoV-2 mutant spectra sampled from infected patients, within the framework of quasispecies dynamics. The analyses suggest that the information provided by myriads of genomic sequences within infected individuals may have a predictive value of the genomic sequences that acquire epidemiological relevance. Possibilities to reconcile the presence of broad mutant spectra in the large RNA coronavirus genome with its encoding a 3' to 5' exonuclease proofreading-repair activity are considered. Indeterminations in the behavior of individual viral genomes provide a benefit for the survival of the ensemble. We propose that this concept falls in the domain of "stochastic thinking," a notion that applies also to cellular processes, as a means for biological systems to face unexpected needs.

In-Host HEV Quasispecies Evolution Shows the Limits of Mutagenic Antiviral Treatments.

Here, we report the in-host hepatitis E virus (HEV) quasispecies evolution in a chronically infected patient who was treated with three different regimens of ribavirin (RBV) for nearly 6 years. Sequential plasma samples were collected at different time points and subjected to RNA extraction and deep sequencing using the MiSeq Illumina platforms. Specifically, we RT-PCR amplified a single amplicon from the core region located in the open-reading frame 2 (ORF2). At the nucleotide level (genotype), our analysis showed an increase in the number of rare haplotypes and a drastic reduction in the frequency of the master (most represented) sequence during the period when the virus was found to be insensitive to RBV treatment. Contrarily, at the amino acid level (phenotype), our study revealed conservation of the amino acids, which is represented by a high prevalence of the master sequence. Our findings suggest that using mutagenic antivirals concomitant with high viral loads can lead to the selection and proliferation of a rich set of synonymous haplotypes that express the same phenotype. This can also lead to the selection and proliferation of conservative substitutions that express fitness-enhanced phenotypes. These results have important clinical implications, as they suggest that using mutagenic agents as a monotherapy treatment regimen in the absence of sufficiently effective viral inhibitors can result in diversification and proliferation of a highly diverse quasispecies resistant to further treatment. Therefore, such approaches should be avoided whenever possible.

The influence of immune cells on the existence of virus quasi-species.

This article investigate a nonlocal reaction-diffusion system of equations modeling virus distribution with respect to their genotypes in the interaction with the immune response. This study demonstrates the existence of pulse solutions corresponding to virus quasi-species. The proof is based on the Leray-Schauder method, which relies on the topological degree for elliptic operators in unbounded domains and a priori estimates of solutions. Furthermore, linear stability analysis of a spatially homogeneous stationary solution identifies the critical conditions for the emergence of spatial and spatiotemporal structures. Finally, numerical simulations are used to illustrate nonlinear dynamics and pattern formation in the nonlocal model.

From viral democratic genomes to viral wild bunch of quasispecies.

The tremendous majority of RNA genomes from pathogenic viruses analyzed and deposited in databases are consensus or "democratic" genomes. They represent the genomes most frequently found in the clinical samples of patients but do not account for the huge genetic diversity of coexisting genomes, which is better described as quasispecies. A viral quasispecies is defined as the dynamic distribution of nonidentical but closely related mutants, variants, recombinant, or reassortant viral genomes. Viral quasispecies have collective behavior and dynamics and are the subject of internal interactions that comprise interference, complementation, or cooperation. In the setting of SARS-CoV-2 infection, intrahost SARS-CoV-2 genetic diversity was recently notably reported for immunocompromised, chronically infected patients, for patients treated with monoclonal antibodies targeting the viral spike protein, and for different body compartments of a single patient. A question that deserves attention is whether such diversity is generated postinfection from a clonal genome in response to selection pressure or is already present at the time of infection as a quasispecies. In the present review, we summarize the data supporting that hosts are infected by a "wild bunch" of viruses rather than by multiple virions sharing the same genome. Each virion in the "wild bunch" may have different virulence and tissue tropisms. As the number of viruses replicated during host infections is huge, a viral quasispecies at any time of infection is wide and is also influenced by host-specific selection pressure after infection, which accounts for the difficulty in deciphering and predicting the appearance of more fit variants and the evolution of epidemics of novel RNA viruses.

Genetic diversity of hepatitis B virus quasispecies in different biological compartments reveals distinct genotypes.

The selection pressure imposed by the host immune system impacts hepatitis B virus (HBV) quasispecies variability. This study evaluates HBV genetic diversity in different biological fluids. Twenty paired serum, oral fluid, and DBS samples from chronic HBV carriers were analyzed using both Sanger and next generation sequencing (NGS). The mean HBV viral load in serum was 5.19 ± 4.3 log IU/mL (median 5.29, IQR 3.01-7.93). Genotype distribution was: HBV/A1 55% (11/20), A2 15% (3/20), D3 10% (2/20), F2 15% (3/20), and F4 5% (1/20). Genotype agreement between serum and oral fluid was 100% (genetic distances 0.0-0.006), while that between serum and DBS was 80% (genetic distances 0.0-0.115). Two individuals presented discordant genotypes in serum and DBS. Minor population analysis revealed a mixed population. All samples displayed mutations in polymerase and/or surface genes. Major population analysis of the polymerase pointed to positions H122 and M129 as the most polymorphic (≥ 75% variability), followed by V163 (55%) and I253 (50%). Neither Sanger nor NGS detected any antiviral primary resistance mutations in the major populations. Minor population analysis, however, demonstrated the rtM204I resistance mutation in all individuals, ranging from 2.8 to 7.5% in serum, 2.5 to 6.3% in oral fluid, and 3.6 to 7.2% in DBS. This study demonstrated that different fluids can be used to assess HBV diversity, nonetheless, genotypic differences according to biological compartments can be observed.

Long-Read Sequencing with Hierarchical Clustering for Antiretroviral Resistance Profiling of Mixed Human Immunodeficiency Virus Quasispecies.

BACKGROUND: HIV infections often develop drug resistance mutations (DRMs), which can increase the risk of virological failure. However, it has been difficult to determine if minor mutations occur in the same genome or in different virions using Sanger sequencing and short-read sequencing methods. Oxford Nanopore Technologies (ONT) sequencing may improve antiretroviral resistance profiling by allowing for long-read clustering. METHODS: A new ONT sequencing-based method for profiling DRMs in HIV quasispecies was developed and validated. The method used hierarchical clustering of long amplicons that cover regions associated with different types of antiretroviral drugs. A gradient series of an HIV plasmid and 2 plasma samples was prepared to validate the clustering performance. The ONT results were compared to those obtained with Sanger sequencing and Illumina sequencing in 77 HIV-positive plasma samples to evaluate the diagnostic performance. RESULTS: In the validation study, the abundance of detected quasispecies was concordant with the predicted result with the R2 of > 0.99. During the diagnostic evaluation, 59/77 samples were successfully sequenced for DRMs. Among 18 failed samples, 17 were below the limit of detection of 303.9 copies/µL. Based on the receiver operating characteristic analysis, the ONT workflow achieved an F1 score of 0.96 with a cutoff of 0.4 variant allele frequency. Four cases were found to have quasispecies with DRMs, in which 2 harbored quasispecies with more than one class of DRMs. Treatment modifications were recommended for these cases. CONCLUSIONS: Long-read sequencing coupled with hierarchical clustering could differentiate the quasispecies resistance profiles in HIV-infected samples, providing a clearer picture for medical care.

Ongoing HIV-1 evolution and reservoir reseeding in two elite controllers with genetically diverse peripheral proviral quasispecies.

BACKGROUND: Elite controllers (EC) are human immunodeficiency virus (HIV)-positive individuals who can maintain low viral loads for extended periods without antiretroviral therapy due to multifactorial and individual characteristics. Most have a small HIV-1 reservoir composed of identical proviral sequences maintained by clonal expansion of infected CD4+ T cells. However, some have a more diverse peripheral blood mononuclear cell (PBMC)-associated HIV-1 reservoir with unique sequences. OBJECTIVES: To understand the turnover dynamics of the PBMC-associated viral quasispecies in ECs with relatively diverse circulating proviral reservoirs. METHODS: We performed single genome amplification of the env gene at three time points during six years in two EC with high intra-host HIV DNA diversity. FINDINGS: Both EC displayed quite diverse PBMCs-associated viral quasispecies (mean env diversity = 1.9-4.1%) across all time-points comprising both identical proviruses that are probably clonally expanded and unique proviruses with evidence of ongoing evolution. HIV-1 env glycosylation pattern suggests that ancestral and evolving proviruses may display different phenotypes of resistance to broadly neutralising antibodies consistent with persistent immune pressure. Evolving viruses may progressively replace the ancestral ones or may remain as minor variants in the circulating proviral population. MAIN CONCLUSIONS: These findings support that the high intra-host HIV-1 diversity of some EC resulted from long-term persistence of archival proviruses combined with the continuous reservoir's reseeding and low, but measurable, viral evolution despite undetectable viremia.

Traceability of SARS-CoV-2 transmission through quasispecies analysis.

During COVID-19 pandemic, consensus genomic sequences were used for rapidly monitor the spread of the virus worldwide. However, less attention was paid to intrahost genetic diversity. In fact, in the infected host, SARS-CoV-2 consists in an ensemble of replicating and closely related viral variants so-called quasispecies. Here we show that intrahost single nucleotide variants (iSNVs) represent a target for contact tracing analysis. Our data indicate that in the acute phase of infection, in highly likely transmission links, the number of viral particles transmitted from one host to another (bottleneck size) is large enough to propagate iSNVs among individuals. Furthermore, we demonstrate that, during SARS-CoV-2 outbreaks when the consensus sequences are identical, it is possible to reconstruct the transmission chains by genomic investigations of iSNVs. Specifically, we found that it is possible to identify transmission chains by limiting the analysis of iSNVs to only three well-conserved genes, namely nsp2, ORF3, and ORF7.

Nucleos(t)ide analogues altered quasispecies composition of hepatitis B virus (HBV)-resistant mutations in serum HBV DNA and serum HBV RNA.

Serum hepatitis B virus (HBV) RNA is a new serological indicator reflecting viral replication with good clinical application prospects. This study aimed to clarify the dynamic changes of serum HBV RNA levels and the quasispecies of HBV RNA virus-like particles in nucleos(t)ide analogues (NAs)-experienced chronic hepatitis B (CHB) patients harboring NAs-resistant mutations and their identifiable effects on NAs resistance. We included CHB patients who were on long-term NAs treatment and with HBV DNA rebound. The longitudinally dynamics of serum HBV RNA levels were quantitatively detected, and the quasispecies differences between serum HBV DNA and serum HBV RNA were compared by high-throughput sequencing. The effect of NAs concentration pressure on altering the resistance mutations quasispecies proportion of HBV DNA and HBV RNA in cell supernatant was analyzed in vitro. A total of 447 serum samples from 36 CHB patients treated with NAs were collected. The median follow-up period was 47 months (about 4 years), and the longest follow-up period was 117 months (about 10 years). Our results showed that HBV RNA could reflect virological breakthrough in 23 (64%, 23/36) patients, and serum HBV RNA rebound earlier than HBV DNA in 12 (52%, 12/23) patients. However, serum HBV RNA remained at a consistently high level and did not fluctuate significantly with the HBV DNA rebound in 6 of 36 patients. In addition, serum HBV RNA was not consistently detectable in 7 of the 36 patients, and their serum HBV RNA was undetectable even after HBV DNA had rebounded. The proportion of drug-resistant mutations in HBV DNA was higher than that of HBV RNA by high-throughput sequencing. The results of in vitro experiments showed that the viral strains with drug-resistant mutation in HBV DNA in cell supernatants gradually become the dominant strains with the increase of NAs concentrations. Serum HBV RNA levels can reflect virological breakthrough in most NAs- treated CHB patients, but there are certain limitations. NAs alter the quasispecies composition of serum HBV DNA and serum HBV RNA, resulting in a higher detection rate of drug-resistant mutations in serum HBV DNA than in serum HBV RNA.

Quantifying In-Host Quasispecies Evolution.

What takes decades, centuries or millennia to happen with a natural ecosystem, it takes only days, weeks or months with a replicating viral quasispecies in a host, especially when under treatment. Some methods to quantify the evolution of a quasispecies are introduced and discussed, along with simple simulated examples to help in the interpretation and understanding of the results. The proposed methods treat the molecules in a quasispecies as individuals of competing species in an ecosystem, where the haplotypes are the competing species, and the ecosystem is the quasispecies in a host, and the evolution of the system is quantified by monitoring changes in haplotype frequencies. The correlation between the proposed indices is also discussed, and the R code used to generate the simulations, the data and the plots is provided. The virtues of the proposed indices are finally shown on a clinical case.

Cellular senescence as a source of SARS-CoV-2 quasispecies.

In-depth analysis of SARS-CoV-2 biology and pathogenesis is rapidly unraveling the mechanisms through which the virus induces all aspects of COVID-19 pathology. Emergence of hundreds of variants and several important variants of concern has focused research on the mechanistic elucidation of virus mutagenesis. RNA viruses evolve quickly either through the error-prone polymerase or the RNA-editing machinery of the cell. In this review, we are discussing the links between cellular senescence, a natural aging process that has been recently linked to SARS-CoV-2 infection, and virus mutagenesis through the RNA-editing enzymes APOBEC. The action of APOBEC, enhanced by cellular senescence, is hypothesized to assist the emergence of novel variants, called quasispecies, within a cell or organism. These variants when introduced to the community may lead to the generation of a variant of concern, depending on fitness and transmissibility of the new genome. Such a mechanism of virus evolution may highlight the importance of inhibitors of cellular senescence during SARS-CoV-2 clinical treatment.